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We present an analytical expression for the response of a transient spectrum to a single- 
Cooper-pair box biased by a classical voltage and irradiated by a single-mode quantized 
field. The exact solution of the model is obtained, by means of which we analyze the an- 



alytic form of the fluorescence spectrum using the transitions among the dressed states of 
the system. An interesting relation between the fluorescence spectrum and the dynamical 
evolution is found when the initial field states are prepared in binomial states. 



> 1 Introduction 

o 

For quantum information science and technologies, it is crucial to build the fundamental 
quantum logic gates pQ. Together with the basic single bit logic gates, the non-trivial 
two bit gates constitutes the fundamental blocks for the quantum network of quantum 
computing. The present lack of a current standard based on quantum devices has inspired 
several attempts to manipulate single electrons, or Cooper pairs, where the rate of particle 
transfer is controlled by an external frequency. Various superconducting nanocircuits have 
qh been proposed as quantum bits (qubits) for a quantum computer In principle, any 

two-state quantum system works as a qubit, the fundamental unit of quantum information. 
^ ■ However, only a few real physical systems have worked as qubits, because of requirements 

of a long coherent time and operability. Among various physical realizations, such as 
ions traps, QED cavities, quantum dots and NMR etc., superconductors with Josephson 
junctions offer one of the most promising platforms for realizing quantum computation 

CO □ E3 

The single-Cooper-pair transistor pHE], is composed of two ultrasmall tunnel junc- 
tions in series forming an island. Transport through the single-Cooper-pair transistor 
depends on the electrostatic energy required to charge the island, as in the single-electron 
transistor, and also on the Josephson coupling across the junctions. Cooper-pair boxes 
are one of the prominent candidates for qubits in a quantum computer. Recent experi- 
ments jH] have revealed quantum coherent oscillations in two CPBs coupled capacitively 
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and demonstrated the feasibility of a conditional gate as well as creating macroscopic 
entangled states. Scalable quantum-computing schemes |Hj have been proposed based on 
charge qubits. In architectures based on Josephson junctions coupled to resonators, the 
resonators store single qubit states, transfer states from one Josephson junction to an- 
other, entangle two or more Josephson junctions, and mediate two-qubit quantum logic. 
In effect, the resonators are the quantum computational analog of the classical memory 
and bus elements. 

In this paper we deal with the problem of the interaction between a single-mode 
quantized field and a single-Cooper-pair box biased by a classical voltage. Despite the 
complexity of the problem, we obtain a quite simple master equation that is valid for 
arbitrary values of the Rabi frequency and the detuning. We apply the Fourier transform 
of the time averaged dipole-dipole correlation function to calculate the fluorescence spec- 
trum, assuming that the electromagnetic field is initially in a binomial state. We find that 
the detuning changes considerably the shape of the resonance fluorescence spectrum and 
leads to novel spectral features. The organization of this paper is as follows: in section 2 
we introduce the model and give exact expression for the unitary operator. In section 3 we 
employ the analytical results obtained in section 2 and by using the finite double- Fourier 
transform of the two-time field correlation function we find an analytical expression for 
the spectrum. Finally, we summarize the results in section 4. 



Several schemes have been proposed for implementing quantum computer hardware in 
solid state quantum electronics. These schemes use electric charge, magnetic flux, su- 
perconducting phase, electron spin, or nuclear spin as the information bearing degree of 
freedom [IT]. In this paper, we consider an example of a realistic system, fabricated by 
the present day technology. We consider a superconducting box with a low-capacitance 
Josephson junction (with the capacitance Cj and Josephson energy Ej), biased by a clas- 
sical voltage source V g through a gate capacitance C g and placed inside a single-mode 
microwave cavity. Suppose the gate capacitance is screened from the quantized radiation 
field, then the junction-field Hamiltonian, in the interaction picture, can be written as 



2 The model 



H23 



(Q - C g V g - CjVf 



Ej cos (f>. 
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where the relevant conjugate variables are the charge Q = 2Ne on the island (where N is 
the number of Cooper-pairs) and the phase difference <fi across the junction. The radiation 
field is to produce an alternating electric field of the same frequency across the junction, 
and V is the effective voltage difference produced by the microwave across the junction. 
We assume that the dimension of the device is much smaller than the wavelength of the 
applied quantized microwave (which is a realistic assumption), so the spatial variation in 
the electric field is negligible. We also assume that the field is linearly polarized, and is 
taken perpendicular to the plane of electrodes, then V may be written down as 



where $ and are the creation and annihilation operators of the microwave with fre- 
quency 10. Cf is the capacitance parameter, which depends on the thickness of the 
junction, the relative dielectric constant of the thin insulating barrier, and the dimension 
of the cavity. In this paper we consider the case where the charging energy with scale 

e 2 

E C = -J (Cg + Cj) , 

dominates over the Josephson coupling energy E T , and concentrate on the value Vg = ^f- , 
so that only the low-energy charge states N = and N — 1 are relevant. In this case the 
Hamiltonian in a basis of the charge state ||) and ||) reduces to a two-state form. In a 
spin- 1/2 language [T3] 



H in = E c [l+ -i-V 2 - -Eja x + 2E c -^Va z . (3) 



- \ejo x + 2E C ^ 
e 2 J 2 e 

where a x and a z are the Pauli matrices in the pseudo-spin basis. 

The time development of the state vector \ip{t)) of the system is postulated to be 
determined by Schrodinger equation 

ihj t m)) = H m \m)- (4) 

The solution of equation (4) can be written as \4>(t)) = U t |^(0)), where U t is a time evolu- 
tion operator which depends on the Hamiltonian of the system. Since \ip{0)) is arbitrary, 
the time evolution operator obeys ihdU t /dt = H in U t . Integrating this equation, gives U t 
= exp J Hindtj . If the system is conservative and H{ n is explicitly independent of 

time, then U t reduces to U t = exp {^H in t 
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One may, also, assume that the two eigenstates l^^) are known, along with their 

(n) 

corresponding eigenenergies T^.Then we can write the time evolution operator as 

oo 

U t = J2 {exp(-2T ( , n) t)|$S n) )($S n) | +exp(-^TL n) t)|^ n) )(^ n) |} . (5) 

n=0 

In order to find explicit forms of the eigenvalues and eigenvectors, we consider the weak 
quantized radiation field and then we can neglect the term containing V 2 in equation (J3J) 
and using the rotating wave approximations we obtain 

Ttf> = o(n+l/2)±^, 



We denote by A = Ej — uj the detuning between the Josephson energy and cavity field 
frequency. | ^*X2 ) are given by 



1 1 eCj^TT) V u> 



KT) 



, -i ( (Cj + C g ) (-A + 2/i w ) [2KCj x , 
+ sin (tan [ ] ] |„+ l. |,. (7) 



1 1 eCj^TY) V u 



-i({Cj + c g )(-A + 2 fin ) Fine j 



kt) 



We devote the next section to investigate the general structure of the dipole-dipole cor- 
relation function in terms of the eigenstates and eigenvalues of the system, and evaluate 
the fluorescence spectrum for input binomial states. 



3 The transient spectrum 

For the calculation of the spectrum, we consider the time evolution of the off-diagonal 
density matrix elements of the field while the diagonal density matrix elements remain 
stationary. In this section we derive the physical transient spectrum S{y) by calculating 
the Fourier transform of the time averaged dipole-dipole correlation function (cr+(t + 
r)er_(i)), weighted by the detector response function {ijj{jS)\a + {t + r)cr_(i) |^>(0)) where 
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\ip{0)) is the initial state of the considered system. Then, the transient spectrum is given 
by the expression [T3j 



S(u) = Re\ drexp[-ii/r-7r](V(0)|<T + (t + r)<7_(*)|V'(0))dr , (9) 



where 7 is the detector width and v is the spectrum frequency. 

Using equations (5)-(8), the time evolution of the states can be expressed as 

U t \n,e) = A(n,t)\n,e) + B(n,t)\n + k,g), 

U t \n,g) = A(n - k,t)\n,g) + B(n - k,t)\n - k,e), (10) 

where 



A(n , t) = sin * 1 tair > ( (gLi^W/O^ ) ) exp (_„<»> 



B(n,t) = Isin2(tan-'( Cj + C »)(- A + 2 ''»)-/^ 



2 \ \ eCj^TT; 



X 



exp ( -itTf ) ) - exp (-ztTl n) ) . (12) 



By using the above equations, the correlation function (er+(£ + t)<t_(£)) can be evaluated 
in terms of the coefficients A(n, t) and B(n, t), as follows (the box initially prepared in the 

M 

excited state i.e. the initial state of the system is assumed to be |^(0)) = ^2 P n (v) \ n i e ) 1 

71=0 

where (3 n (if) will be defined latter) 

00 

(a+(t + T)&-(t)) = J2(Pn) 2 (n,e\a + (t + T)a_(t)\n,e) 

n=0 

00 

= J2(f 3 n) 2 Mn,t)A(n-k,r)A*(n,t-r). (13) 

n=0 

The Fourier transform of the time averaged dipole-dipole correlation, which is directly 
related to the fluorescence spectrum with the identification of 7 as the width associated 
with the detector After carrying out the various operations we get 



7 2 + (z/ _ T +)2 7 2 + (i ,_ T -)2 
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+EW)'[ 



7 sin 4 0„ cos 2 



n—k 



n=l 



7 2 + (^ + T+_ fc -T 



+ 



7 sin 4 0„ sin 2 6> n _ fc 



7 cos 4 # n cos 2 6* 



n—k 



+ 



n/ 7 2 + (^ + T n _ fc 
7 cos 4 9 n sin 2 6> n _ fc , 



n I 



7 2 + (i/ + T+_ fc - T-)2 7 2 + + T'_ fe - T") 



(14) 



where (/3^) 2 is the initial photon number distribution and n is 



0, 



tan 



1; 



2HC, 
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Thus the time averaged spectrum consists of resonant structures which arise from transi- 
tions among different dressed states. The final structure of the time averaged spectrum 
will depend on the form of the input photon distribution. As the cavity field starts to 
interact with the Cooper pair the initial photon number distribution starts to change. 
Due to the quantum interference between component states the oscillations in the cavity 
field become to be composed of two component states. The situation that has just been 
described is depicted in figures 1-5. 

1, 
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Figure 1: The evolution of the function S(v) in a perfect cavity as a function of (v — uj)/X 
for an input binomial state, with the r] = 0, M — > oo(i.e.r]M = \a\ 2 where \a\ 2 = 10 (solid 
line) and \a\ = 1( dotted line). Calculations assume detector width 7 = O.lg, and the 
detuning parameter A has zero value. 

In figure 1 we discuss the spectrum for the situation of a coherent field with A = 0. 
We show the evolution of this spectrum S(v) as a function of (v — u)/X where A = 



6 



eC.i 



2{c,+c ) y 2^7- We cons ider the initial state of the field is a binomial state. The binomial 
states which has been introduced by Stoler, Saleh and Teich in [2], interpolate between the 
most nonclassical states, such as number states and coherent states, and reduce to them 
in two different limits. Some of their properties [2, 3, 4], methods of generation [2, 3, 5], 
as well as their interaction with atoms [6], have been investigated in the literature. The 
binomial state is defined as a linear superposition of number states in an M— dimensional 
subspace 

M 

\v,M) = J2Pn(v)\n), (15) 

n=0 

where the state |n) is the number state of the field mode, 77 is a complex number with the 
absolute value between and 1, M is a positive integer, and 




{rf 1 + (1 — rj) n ~ M ). (16) 



The name binomial state comes from the fact that their photon distribution is simply a 
binomial distribution with probability 77. The binomial state is a linear combination of 
M + 1 number states with coefficients chosen such that the photon-counting probability 
distribution is binomial with mean photon M. This state can produce, under certain 
choices of the parameters rj and M, the number state |M), the vacuum state |0) and the 
coherent state. When rj — > 0, M — > 00 with r/M = a 2 fixed (a real constant), \rj : M) 
reduces to the coherent states which correspond to the Poisson distribution. Here we 
consider 77 = 0, M — > 00 and \a\ 2 = 10 (solid line) and |a| = 1 (dotted line). Calculations 
assume that the detector width 7 = 0.1 A, and the detuning parameter A has zero value. 
The slit is adjusted so that the light is collected from Cooper-pair box which have been in 
the field for times ranging from up to 7. From an initially broad featureless spectrum, 
the central peak and the sidebands emerge rather quickly. As the observation region is 
lengthened, all the components get taller and narrower. Then for about (v — uo) ps ±9A, 
the transient spectrum has narrowed enough and tends to zero. In the limit of a very large 
mean photon number and at exactly resonant field, it is possible to simplify the expression 
for the spectrum S(v) and illustrate explicitly the manner in which the sidebands and the 
central peaks narrow with increasing (v — uo). 

It would be of interest to pay attention to the physical transient spectrum for binomial 
states. In order to do that, we use different values of both r] and M (say M = 3 and 
r) = 0.7, 0.1). This basis is overcomplete and many states in the Hilbert space can 
be expanded in it. The behaviors of S(u) is notably different from those observed in 
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Figure 2: The evolution of the function S(v) in a perfect cavity as a function of (v — uj)/\ 
for an input binomial state, with M = 3, rj = 0.7 for (solid line) and rj = 0.1 for (dotted 
line). Calculations assume that, the detector width 7 = O.lg, and the detuning parameter 
A has zero value. 

the coherent state case. Indeed, in the fully connected system considered here, S(u) 
is extremum at ±1(±0.4) when rj = 0.7(0.1) whereas in the coherent state case, figure 
1, S(v) is extremum at « 0. In addition, the scaling behavior of the spectrum and of 
its derivative are different in both cases. Therefore, when speaking about the physical 
transient spectrum sensitivity achieved with different frequencies, it is necessary to specify 
the kind of initial state of the field. 

Our aim is now to evaluate the response of the physical transient spectrum due to 
the mean photon M. It is interesting to note that in a large value of M, we observe 
a similar behavior to that obtained in the coherent state, see figure 3. Now, we would 
like to shed some light on the spectrum behavior when the detuning differs from zero. 
The transient spectrum under these conditions can be asymmetric and ultimately the 
central component and one of the Rabi sidebands can vanish despite the fact that the 
quadrature-noise spectrum exhibits a significant amount of noise at these frequencies. 
The asymmetry arises from the stimulated emission induced between the dressed states 
by the binomial field. 

These asymmetries consist of an enhanced sideband on the atomic resonance frequency 
side of the central peak and more pronounced oscillations between the central peak and 
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Figure 3: The same as in figure 2 but M = 30. 
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Figure 4: The same as in figure 2 but A = g. 



the enhanced sideband than on the opposite side of the central peak. Also, we observe a 
slight displacement in the location of the central maximum from the applied field frequency 
toward the atomic frequency. One can see that the large value of the detuning parameter 
gives a disappearance of one peak and large displacement in the location of the maximum 
value of the spectrum. 
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Figure 5: The same as in figure 2 but A = 2g. 

We can conclude that the effect of the detuning on the spectrum of the emitted light 
is twofold. The first effect is the shift of the spectrum to the left or to the right depending 
on the sign of A. The second effect is the dependence of the amplitude of the peaks on 
A. This dependence leads to the fact that, in the far off-resonance limit, only one peak 
survives. 

Finally we may point out that the present lack of a current standard based on quantum 
devices has inspired several attempts to manipulate single electrons, or Cooper pairs, 
where the rate of particle transfer is controlled by an external frequency [16| . In the 
experiments described in [T7] a long array of Josephson junctions with an external signal 
applied to a gate, which is capacitively coupled to the middle of the array has been 
used. Theoretically many systems can act as a qubit, but the realization of it is difficult. 
Nakamura et al. ^H] have shown in their experiments with single Cooper pair box, that 
in a kind of metallic island structures the oscillations between eigenstates of the system 
last at least a few nanoseconds. The states are characterized by the number of Cooper 
pairs in the box and (quantum) manipulation of this number is of basic importance for 
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the production of viable qubits. This might be a very important issue when thinking of 
the limitations on the preparation and read-out of the states. 

4 Conclusion 

We have analyzed the physical transient spectrum of a single Cooper-pair box, which is 
biased by a classical voltage and irradiated by a single-mode quantized field. We emphasize 
the fact that the proper expression for the emission spectrum which is derived in this 
paper can be measured in a realistic experiment. This spectrum has been obtained not 
only as a function of the atomic and field parameters, but also as a function of 7 which is 
available to the experimenter (7 the width associated with the detector). Thus it correctly 
incorporates the possible effects of a finite observation interval, observations made close 
to the point where the interaction was turned on, and arbitrary initial conditions for 
the Cooper pair box at the start. Among the reasons for the interest in considering 
nonclassical effects of the binomial states, we may mention the following reasons: the 
experimental work shows that nonclassical effects serve as a test of the quantum nature of 
light and nonclassical behavior of light is usually connected with a noise reduction below 
a standard limit (e.g. the shotnoise limit). In particular, we have explored the influence of 
various parameters of the system on the emission spectrum of the output field statistics. 
We have used the finite double-Fourier transform of the two-time field correlation function 
to find an analytical expression for the spectrum. The spectrum in the cavity for the 
initial binomial states is studied. Such systems are potentially interesting for their ability 
to process information in a novel way and might find application in models of quantum 
logic gates. The phenomenon of oscillations in the field spectrum has been shown. It 
is observed that the symmetry shown in the resonant case for the spectra is no longer 
present once the detuning is added. 
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